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The decay channel ip' ■n'^ -k' J /ij){J /ip — > 7Pp) is studied using a sample of 1.06 x 10® ip' events 
collected by the BESIII experiment at BEPCII. A strong enhancement at threshold is observed in 
the pp invariant mass spectrum. The enhancement can be fit with an 5*- wave Breit-Wigner resonance 
function with a resulting peak mass of M = 186lli3 (stat)+2g (syst) MeV/c^ and a narrow width 
that is r < 38 MeV/c^ at the 90% confidence level. These results are consistent with published 
BESIl results. These mass and width values do not match with those of any known meson resonance. 

PACS numbers: 13.85.Hd, 25.75.Gz 



An anomalously strong pp mass threshold enhance- 
ment was observed by the BESII experiment in the 
radiative decay process J/V' jpp [H- In ref. [1] 
it was noted that when an S'-wave Breit-Wigner res- 
onance function is fitted to the pp mass distribution, 
the peak mass is below the pp mass threshold at M = 
1859tio (stat)+25 (syst) MeV/c^ and the total width is 
r < 30 MeV/c^ (at the 90% C.L.). An interesting fea- 
ture of this enhancement is that corresponding structures 
are not observed in near-threshold »» cross section mea- 
surements, in B- meson decays @, 0, in radiative tp' or 
T — > 'fpp decays [1, or in J/tp — ojpp decays Q. 
These non-observations disfavor the attribution of the 
mass-threshold enhancement, which is uniquely observed 
in the J/ip — jpp decay process, to the pure effects of 
pp final state interactions (FSI). 

This experimental observation stimulated a number of 
theoretical speculations 0, H, S 113, EH 113 ■ One of these 
is the intriguing suggestion that it is an example of a 
pp bound state, sometimes called baryonium [1^, which 
has a long history and has been the subject of many 
experimental searches [ijj. 



In this letter we report a study of the pp mass spec- 
trum in the threshold region in the decay process ip' — > 
TT^ir^ J/ip(J/ip — >■ jpp). The analysis uses a sample of 
1.06 X 10^ V' events accumulated by the upgraded Beijing 
Spectrometer (BESIII) located at the Beijing Electron- 
Positron Colhder (BEPCII) at the Beijing Institute of 
High Energy Physics. 

BEPCII/BESIII |15|] is a major upgrade of the BE- 
SII experiment at the BEPC accelerator [l^. The de- 
sign peak luminosity of the double-ring e^e~ collider, 
BEPCII, is 10^3 cm-^s-i at a beam current of 0.93 A. 
The BESIII detector with a geometrical acceptance of 
93% of 47r, consists of the following main components: 1) 
a small-celled, helium-based main draft chamber (MDC) 
with 43 layers. The average single wire resolution is 135 
fim, and the momentum resolution for 1 GeV charged 
particles in a 1 T magnetic field is 0.5%; 2) an electromag- 
netic calorimeter (EMC) made of 6240 Csl (Tl) crystals 
arranged in a cylindrical shape (barrel) plus two endcaps. 
For 1.0 GeV photons, the energy resolution is 2.5% in the 
barrel and 5% in the endcaps, and the position resolu- 
tion is 6 mm in the barrel and 9 mm in the endcaps; 3) a 
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Time-Of-Flight system (TOF) for particle identification 
composed of a barrel part made of two layers with 88 
pieces of 5 cm thick, 2.4 m long plastic scintillators in 
each layer, and two endcaps with 96 fan-shaped, 5 cm 
thick, plastic scintillators in each endcap. The time res- 
olution is 80 ps in the barrel, and 110 ps in the endcaps, 
corresponding to a 2cr K/tt separation for momenta up to 
about 1.0 GcV; 4) a muon chamber system (MUC) made 
of 1000 m^ of Resistive Plate Chambers (RFC) arranged 
in 9 layers in the barrel and 8 layers in the endcaps and 
incorporated in the return iron of the superconducting 
magnet. The position resolution is about 2 cm. 

Candidate ip' — ^ ■k^'k~ J /ipl J /ip ipp) events are 
required to have at least one photon and four charged 
tracks within the polar angle region jcos^j < 0.93 and 
a total net charge of zero. The TOF and dE/dx infor- 
mation are combined to form particle identification con- 
fidence levels for the tt, K, and p hypotheses; each track 
is assigned to the particle type that corresponds to the 
hypothesis with the highest confidence level. Selected 
events are required to have both an identified proton and 
an identified anti-proton and no particle identification is 
required for the two remaining tracks. Candidate pho- 
tons are required to have an energy deposit that is at 
least 25 MeV in the barrel EMC (| cos 6*1 < 0.8)) and 
50 MeV in the endcap EMCs (0.86 < |cos6i| < 0.92), 
and be isolated from the anti-proton track by more than 
30° due to the strong annihilation of anti-protons, and 
from all other charged tracks by more than 10°. EMC 
timing requirements suppress electronic noise and energy 
deposits unrelated to the event. 

Candidate J/^ signals are identified by the invari- 
ant mass recoiling against the tt^tt^' pair, \M^+^- — 

iT^,j/4>\ < 0.006 GeV/c^. Further requirements of 
\Umiss\ < 0.05 GeV, where Umiss = {Emiss - 
< 0.0005 (GeV/c)2 



\Pmiss\), and Pf^ 



where P^^ 



4|Pmiss| sin 6*7/2, are imposed to suppress backgrounds 
from multi-photon events. Here Emiss 

and Pjriiss are, 

respectively, the missing energy and momentum of all 
charged particles, and 9-^ is the angle between the missing 
momentum and the photon direction. The requirement 
\M^+„-pp — TO^'I > 0.03 GeV/c^ is used to reduce the 
background from tp' TT'^TT~pp. 

Events that remain after these selection requirements 
are subjected to a four-constraint energy- momentum 
conservation kinematic fit to the hypothesis ip' — > 
jTT~^TT~pp. For events with more than one 7 candidate, 
the combination with the smallest is chosen. Events 
with < 100 are selected. Since the detection efficien- 
cies for data and Monte Carlo (MC) simulated events 
are consistent for protons and antiprotons with momenta 
p > 0.3 GeV/c, while differences occur for lower momen- 
tum tracks, we reject events with pp < 0.3 GeV/c or 
Pp < 0.3 GeV/c. 

Fig. [ija) shows the pp invariant mass distribution for 
surviving events. The distribution features a peak near 
Mpp — 2.98 GeV/c^ that is consistent in mass, width, 
and yield with expectations for tp' tt+tt^ J/-0( J/-0 
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FIG 

ip' — > ■k'^-k~ J /ip{J /ip — >■ 7pp) candidate events, (a) The pp 
invariant mass spectrum; the open histogram is data and the 
hatched histogram is from a -0' — >■ 'k'^tx~ J /■ip{J /^p — > 7pp) 
phase-space MC events(with arbitrary normalization), (b) 
An A/'^(7p) (horizontal) versus M^ljp) (vertical) Dalitz plot 
for the selected events. 



IVciVc ~^ PP)j a broad enhancement around Mpp ^ 
2.2 GeV/c^, and a prominent low-mass peak at the pp 
mass threshold, similar to that reported by BESII [1]. 
The Dalitz plot for selected events is shown in Fig. [T](b), 
where a band corresponding to the threshold enhance- 
ment is evident in the upper right corner. 

Potential background processes are studied with an in- 
clusive MC sample of 1 x 10* ip' events generated ac- 
cording to the Lund-Charm model 11 Til and the Particle 
Data Group (PDG) decay tables |l8| . None of back- 
ground sources produce an enhancement at the thresh- 
old region of pp invariant-mass spectrum. The domi- 
nant background is from Tp>' — > tt^tt^ J/ip(J/ip — >■ n'^pp) 
events, with asymmetric tt'^ 77 decays where one of 
the photons has most of the n'^ energy. An exclusive 
MC sample of tp' -k^'k~ J /ip{J /ip 7r°pp), generated 
with a uniform phase space distribution, indicates that 
the level of this background in the selected event sample 
with Mpp - 2mp < 0.3 GeV/c^ is 9% of the total. 

To ensure further that the pp threshold enhancement 
is not due to background, each potential background is 
studied with data. Non-J/tp background are studied us- 
ing J/ip mass-sideband events. For these there is no en- 
hancement and their level of contamination of the se- 
lected event sample is about 2%. The dominant back- 
ground channel, ip' — >■ tt'^tt~ J/ip(J/ip i:^pp), is also 
studied with data. In this case, events with four charged 
tracks, including a proton and antiproton and two oppo- 
sitely charged pions, and with two or more photons are 
selected, and subjected to a four-constraint kinematic fit 
to the ip' — > jjTr'^Tr~pp hypothesis. J/^p and tt'^ signals 



< 



are selected by the requirements \M^+^- 

0.006 GeV/c^ and \M-y-y - m^o\ < 0.008 GeV/c^ (±2cr). 
There is no evidence of a narrow and strong enhancement 
near pp mass threshold region. 

The Mpp invariant mass spectrum in the threshold re- 
gion for the selected n'^pp events is shown in Fig. ^a), 
where no threshold enhancement is evident. The dis- 
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FIG. 2: Thepp mass spectrum near threshold for: (a) selected 
ip' tt'^ty' J/iplJ/tp —>■ Tv^pp) events for the same real data 
sample, (b) phase-space MC ^' -> ■K'^n~ J/ip{J/ip 7Pp) 
events that satisfy the 7pp selection criteria. The smooth 
curves are the results of the fit described in the text. 
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FIG. 3: The pp invariant mass spectrum for the ^' — >■ 
tt'^-k" J/ip{J/ip "fpp) after final event selection. The solid 
curve is the fit result; the dashed curve shows the fitted back- 
ground function, and the dash-dotted curve indicates how the 
acceptance varies with pp invariant mass. 



tribution is well described by a function of the form 
fbkgiS) = N{S^/'^+ai6^/'^+a2S^/^), where 6 = Mpp-2mp 
and the shape parameters ai and 02 are determined from 
a fit to selected jpp events for ij;' TT~^Tr~ J/ip{J/ip 
"ipp) phase-space MC sample shown in Fig.[2j6). 

To characterize the pp threshold mass enhancement, 
we fit it with an acceptance weighted Breit-Wigner (BW) 

2L + 1.3 

function of the form BW{M) oc (jy^2_V^)2_)_jy^iir2 ; where 
r is a constant (determined from fit), q is the proton 
momentum in the pp rest-frame, L is the pp orbital an- 
gular momentum, and k is the photon momentum, to- 
gether with the function fbkg{5) with free normalization 
and constants oi and 02 fixed at the ir'^pp phase-space 
MC values {i.e. the curves shown in Fig. [2]) to repre- 
sent the background from mis-reconstructed ir^pp events 
and a possible non-resonant pp phase-space contribution. 
The BW is multiplied by the MC-determined signal ac- 
ceptance that is corrected for MC and data differences of 
the low momentum tt"*" and tt^ tracking efficiencies. The 
tracking efficiencies determined from data are measured 
using samples of tagged protons and antiprotons from the 
process J/tp — >■ pp7r+7r~. 

The result of a fit using L — Q and confined to the 



Mpp — 2mp < 0.3 GeV/c^ mass region is shown in FiglJ] 
The fit returns a signal yield of SlOlgg (stat), a peak 
mass of M = 18611^3 (stat) MeV/c^ and a width of T = 
± 23 MeV/c2. The fit quality is x^/d.o.f. = 42.6/56. 

In the above-described fit, the phase-space contribu- 
tion is treated as an incoherent background under the 
enhancement. Possible fitting biases near threshold are 
investigated using a set of MC samples that combine 
the signal with a uniform incoherent phase-space back- 
ground. In each MC sample, the mass, width, and num- 
ber of signal events are obtained from a fit using the same 
procedure as that done on the data. The averaged dif- 
ferences between the fitted output and input values are 
taken as a systematic uncertainty associated with a pos- 
sible fitting bias. The r.m.s. of each parameter's bias 
measurements is taken as the statistical error. The sys- 
tematic uncertainties found from varying the bin size and 
the fitting range are also included. The total systematic 
error on the mass is I26 MeV/c^. Including the system- 
atic error, the upper limit on the width is F < 38 MeV/c^ 
at a 90% confidence level. 

We also tried to fit the pp mass spectrum using known 
particle resonances to represent the low-mass peak. 
There are two spin-zero resonances listed in the PDG 
tables in this mass region: the 77(1760) with M^(i76o) = 
1756 ± 9 MeV/c2 and r,,(i760) = 96 ± 70 MeV/c^, and 
the 7r(1800) with M^(i8oo) = 1816 ± 14 MeV/c^ and 
r7r(i800) = 208 ± 12 MeV/c^. A fit with fbkg and an 
acceptance- weighted S'-wave BW function with mass and 
width fixed at the PDG values for the 77(1760) produces 

/d.o.f. — 144.8/56; and using the 7r(1800) parameters 
produces /d.o.f. = 161.7/56. 

In summary, an anomalous strong, near-threshold en- 
hancement in the pp invariant mass distribution is ob- 
served in the decay process of 7/;' — > tt+tt^ J/7/;( J/t/; ^■ 
7pp). If it is fitted with an S'-wave Breit- 
Wigner resonance function, the peak mass is M = 
186lli3 (stat)+26 (syst) MeV/c^ and the width is F < 
38 MeV/c^ at the 90% confidence level. These values 
are consistent with the published BESII results [l|. As 
indicated in Ref. [l^, the pp mass threshold enhance- 
ment may also be fitted with a broad structure (F ~ 
100 MeV/c^) multipHed by an FSI factor in Ref. [H]. 
More precise measurement of the shape of the pp mass 
threshold enhancement and more sophisticated fits such 
as including some model dependent FSI factor in the fit 
will be performed with much higher statistics J/t/; data 
sample collected at BESIII. 
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